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Целью исследования является определение кинематических 
параметров сгибания локтя с помощью искусственных 
нейронных сетей (ИНС). В качестве входов ИНС использова- 
лись параметры поверхностной электромиограммы (пЭМГ). 
Выходом ИНС были кинематические параметры движения: 
направление, угловое перемещение и угловая скорость. У 
каждого испытуемого регистрировались сигналы пЭМГ, сни- 
маемые с бицепса при сгибании локтя без нагрузки. Во время 
эксперимента плечевой и локтевой суставы фиксировались 
пассивным экзоскелетом. Формирование вектора признаков 
для нейронной сети выполнялось с помощью методов спек- 
трального и статистического анализов. Статистический анализ 
во временной области включал в себя определение: дисперсии 
амплитудных значений ПЭМГ, среднеарифметического и 
среднеквадратического значений абсолютных амплитуд 
пЭМГ, частоты пересечения нулевой линии сигналом ПЭМГ. 
В частотной области выполнялся спектральный анализ сигна- 
лов ЭМГ методом быстрого преобразования Фурье. Опреде- 
лялись спектр мощности, средняя частота спектра мощности. 
Установлены параметры пЭМГ, использование которых в 
качестве входных сигналов ИНС обеспечивает наименьшую 
ошибку оценки параметров движения. ИНС обучалась мето- 
дом прямого распространения сигнала и обратного распро- 
странения ошибки. Наилучшие результаты определения кине- 
матических параметров получены при использовании в каче- 
стве входов ИНС средней частоты спектра мощности и общей 
интегральной мощности сигнала пЭМГ. Полученные резуль- 
таты могут быть использованы при разработке систем био- 
электрического управления мехатронными устройствами. 


Ключевые слова: пассивный экзоскелет, параметры движения 
локтя, электромиографический сигнал, статистический анализ, 
спектральный анализ, искусственная нейронная сеть, МаЧаБ, 
ЗипаНикК. 
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Введение. Электромиографический (ЭМГ) сигнал широко используется для управления различными электромехани- 
ческими устройствами: протезами [1, 2], экзоскелетами [3] и другими мехатронными системами. Для управления био- 
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электрическими протезами применяют ЭМГ-сигналы от двигательных волокон периферического нерва с помощью 
имплантируемых электродов. В управлении большинства мехатронных устройств используются параметры поверх- 
ностной электромиограммы (пЭМГ) [4]. 

Скелетные мышцы представляют собой совокупность двигательных единиц (ДЕ). При стимулировании 
нейронным сигналом каждая ДЕ сокращается и генерирует биоэлектрический сигнал, который является суммой по- 
тенциалов действия всех вовлеченных в процесс клеток. Сумма потенциалов называется потенциалом действия от- 
дельной двигательной единицы (ПДОДЕ). Пространственно-временное суммирование ПДОДЕ всех активных двига- 
тельных единиц приводит к возрастанию ЭМГ мышцы [5, 6]. 

Очевидно, мехатронными устройствами можно управлять при условии удаления помех из исходного сигнала 
пЭМГ и применении искусственных нейронных сетей (ИНС) [7-9]. Проблема получения стабильного пЭМГ-сигнала 
решается использованием биоусилителей на базе инструментальных усилителей с большим коэффициентом ослабле- 
ния синфазных сигналов (помех) [1, 10]. 

В статье приведены методика и результаты исследований ПЭМГ, снятых с бицепса, результаты статистиче- 
ских и спектральных анализов, выполненных с целью определения минимально достаточного вектора входных сигна- 
лов ИНС, представлена оптимальная архитектура ИНС для определения кинематических параметров движения локтя 
Методика исследований. В исследованиях принимали участия одиннадцать студентов ДГТУ (8 мужчин и 3 женщи- 
ны), не имевшие травм и заболеваний мышц плеча и предплечья. Среднеарифметические и среднеквадратические зна- 
чения антропометрических параметров испытуемых приведены в табл. 1. 

Таблица 1 
Антропометрические параметры испытуемых 





Объем бицепса /СКО, м 


Пол | Длина плеча/ | Длина предплечья/ Вес/СКО, Рост/СКО, 
СКО*, м СКО, м Состояние 





кг см 





Расслабленное | Напряженное 





Муж. | 0,325/0,0132 0,3457/0,0093 0,305/0,0235 0,345/0,0367 | 83,57/11,41 178,5/4,31 





Жен. | 0,302/0,0189 0,3283/0,0175 0,2833/0,0321 | 0,295/0,0346 59,0/10,14 164,33/4,04 























*СКО — среднеквадратическое отклонение. 











Для единообразного начального взаиморасположения локтевого и плечевого суставов во время регистрации 
пЭМГ использовался пассивный экзоскелет [11]. В процессе эксперимента испытуемые совершали угловое движение 
(сгибание и разгибание) правым предплечьем и кистью в сагиттальной плоскости, проходящей через локтевой и пле- 
чевой суставы. Условно предплечье с кистью до конца вытянутого среднего пальца руки будем считать локтем. Угло- 
вые движения локтем (сгибание и разгибание) относительно локтевого сустава совершались без нагрузки (рис. 1). 
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Рис. 1. Векторная модель сгибания локтя [12]: ОП — ось плеча; ОПр — ось предплечья; Р — приведенный вес руки; 
Ет — сила бицепса; 9 — угол поворота локтевого сустава; а — расстояние от локтевого сустава до точки вставки бицепса 


В лучевой кости (радиальная бугристость); ь— расстояние от локтевого сустава до центра тяжести предплечья; 
с — расстояние от локтевого сустава до центра кисти; а— расстояние от локтевого сустава до точки фиксации бицепса 
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Для записи сигналов ПЭМГ на поверхности бицепса по линии мышечных волокон крепились электроды 
(Аз-АвзСГ) диаметром 12 мм. Расстояние между ними составляло 2 см. Электроды фиксировались с помощью манже- 
ты для обеспечения надежного контакта с кожей предплечья [11]. Место установки обрабатывалась в соответствии с 
рекомендациями [13]. Сигнал пЭМГ усиливался в 2300 раз и подвергался фильтрации активными полосовыми про- 
пускающими фильтрами Баттерворта с частотами среза 6, 5 и 318 Гц. 

Сигналы пЭМГ и гироскопа (угловое перемещение локтя) оцифровывались с частотой 1 кГц и передавались 
через систему сбора информации М15В-6212 в ЭВМ. Данная система имеет встроенный 16-разрядный аналого- 
цифровой преобразователь (АЦП). Сигналы пЭМГ обрабатывались с помощью системы МаЧаб — ЗипаГаикК. 

Целью экспериментальных исследований являлось изучение состава вектора входных сигналов для искус- 
ственной нейронной сети (ИНС), посредством которой предполагается распознавать сигналы ЭМГ и определять коли- 
чественные параметры движения локтя. 

Определение характеристик пЭМГ. Поверхностная электромиограмма является нестационарным случайным про- 
цессом и характеризует состояние мышцы на основе электрического потенциала, генерируемого ею при совершении 
какого-либо действия. 

Для анализа пЭМГ-сигнала во временной области вычислялись статистические характеристики, в частотной 
области использовалось преобразование Фурье для определения основных частотных оценок. Статистический и спек- 
тральный анализы выполнялись для наборов данных, формируемых с помощью скользящего временного окна разме- 
ром 200 мс, без перекрытия предыдущего окна. 

В режиме псевдореального времени определялись следующие статистические параметры [13, 14]. 

Среднее абсолютных значений — параметр, являющийся среднеарифметическим значением биопотенциалов 
(БП) пЭМГ в скользящем окне (тАУ): 

ТАЙ = Иж, (р 
где М — число значений БИ пЭМГ в скользящем окне; хх — К-е значение БП пЭМГ. 


Среднеквадратическое отклонение (5ТО) характеризует степень рассеивания БП пЭМГ относительно изо- 
электрической оси: 





1 — 
5ТР = Е м (жк - ХР, (2) 


где Х — среднеарифметическое значение БП. 
Среднеквадратическое значение биопотенциалов пЭМГ (М5) — это среднеквадратическое значение, вычис- 
ленное в пределах окна: 


ВМ = КУ (3) 


Число пересечений сигналом ПЭМГ нулевой линии (7С) — это показатель, связанный с частотой и сложно- 
стью сигнала. Для устранения влияния низкоамплитудных шумов при обработке сигнала вводился порог = = 20 мВ. 
Значение 7С инкрементируется, если 


({хк > Оихк-1 < 0} или {хх < их: > (и [хх — хи | > =. (4) 
Параметр \/ауеогт Гепо, характеризующий форму сигнала, является мерой его сложности и определялся на 
основе отсчетов в исследуемом окне [8] следующим образом: 


И’азегеп = УМ. |Ахк|,где Ахк = хк — Хкф1. (5) 
Дисперсия (уаг) пЭМГ-сигнала является усредненной характеристикой его вариабельности относительно 
среднего значения: 


1 ум 2 
тат = — Уфа хЕ. 6 
м—т24К=1 ХК (6) 
Дисперсия сигнала в скользящем окне является показателем средней мощности сигнала как функции времени. 


Дисперсия определена относительно нулевой линии пЭМГ. 
Развиваемое мышечное усилие оценивалось посредством параметра оз)е! по следующей зависимости: 


[одре = емй-1 109 к а 
Спектральный анализ сигналов ПЭМГ выполнялся методом быстрого дискретного преобразования Фурье: 
М-1 м 
еда мМЕт, (8) 
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где М — количество компонентов разложения; х„, — значения входного вектора; Х. — выходной вектор, состоящий из 
М комплексных амплитуд синусоидальных сигналов, составляющих исходный сигнал пЭМГ; } — мнимая единица; е — 
основание натурального логарифма; 4 — порядковый номер комплексных амплитуд синусоидальных сигналов. 

Так же определялись общая мощность спектра в окне, средняя частота, медиана, частота спектральной со- 
ставляющей, имеющей наибольшую амплитуду [13]. 

Общая интегральная мощность сигнала Пе определялась как 


— [№ Му Фи Рлн а) 
тёе = Граф ыы РИ, (9) 
где /о, № — частотный интервал спектра мощности пЭМГ сигнала; р(Ё) — спектральная плотность мощности; № — 
число дискретных частот; } — 1-е значение дискретной частоты; ру: — 1-е значение спектральной плотности мощности 


сигнала пЭМГ. 

Изменение спектральных характеристик оценивалось также с помощью определения средней частоты спектра 
мощности сигнала МЕ$Р, которой предполагается определять параметры движения локтя: 
Нл Е ГЫ 

Л ря 

Результаты исследований. Для решения задач распознавания фаз движения локтя были спроектированы три нейрон- 
ные сети, которые содержали нейроны первого слоя, скрытый слой нейронов и выходной слой. Такое решение обу- 
словлено нелинейностью и нестационарностью связей вычисляемых статистических и частотных оценок с параметра- 
ми движения локтя. Обучение сетей выполнялось по алгоритмам прямого распространения сигнала и обратного рас- 
пространения ошибки. Определение наилучшего вектора входных сигналов для каждой из сетей является нетривиаль- 
ной задачей. При ее решении сравнивались результаты обучения нескольких вариантов каждой из сетей при различ- 
ных оценках пЭМГ-сигнала на входе. Наилучший вариант ИНС выбирался по результатам сравнения среднеквадрати- 


ческой ошибки (СКО) при распознавании. 


МЕЗР = Ё (10) 


Формирование входного вектора выполнялось из следующих характеристик и оценок пЭМГ-сигнала: МЕЗР, 
те, ГС, о8)е!, КМ5, 5ТР. 

Задачи распознавания фаз и параметров движения формулировались следующим образом. Для первой сети: 
определить наличие/отсутствие движения локтя. Для второй сети: оценить угол поворота (положение) локтя. Для тре- 
тьей сети: определить направление движения локтя. 

При исследовании различных вариантов ИНС для каждой из сетей скрытый слой составлял 5, 10, 15 или 20 
нейронов. Нейроны входного и выходного слоя имели линейные функции активации. Сигмоидальная функция приме- 
нялась на скрытом слое. 

Результаты моделирования процессов распознавания для трех вариантов нейронных сетей приведены в 
табл. 2, 3, 4. 

В табл. 2 представлены значения среднеквадратической ошибки распознавания наличия/отсутствия движения 
для различных наборов входных значений. Как видно из таблицы, наименьшей ошибкой обладает ИНС с входным 
вектором МЕ$Р, Ге. Выходное значение ИНС, равное 0, соответствует отсутствию движения, равное | — наличию 
движения. 


Таблица 2 
Значение СКО распознавания наличия/отсутствия движения предплечья 





МЕЪЗР, | 2С, | МЕР, |МЕЪЗР, ше, | МЕЗР, [ше, | МЕ$Р, ше, ГС,| МЕЪЗР, тие, СС, 
ие |108)е! | ше, 7С| 10о8вреё | 27С,1ов)е | 10оз)ер, КМ$ |1ов)еь КМ5, 5ТЬ 


Среднеквадратическое отклонение 


Кол-во нейронов 





в скрытом слое 









































э 0,0284 | 0,0313 | 0,0284 0,0299 0,0287 0,0293 0,0286 
10 0,0294 | 0,0304 | 0,0285 0,0298 0,0285 0,0295 0,0291 
15 0,0295 | 0,0291 | 0,0288 0,0290 0,0289 0,0286 0,0288 
20 0,0293 | 0,0313 | 0,0297 0,0328 0,0293 0,0285 0,0298 





Необходимо отметить, что для обучения и проверки качества распознавания использовались различные набо- 
ры данных. Величина ошибки 0,0284 является малой (см. табл. 2), что позволяет говорить о надежном распознавании 
ИНС наличия или отсутствия движения на основе пЭМГ-сигнала (рис. 2). 
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Рис. 2. Выходные значения ИНС при определении наличия движения локтя: 
соотношение между выходом ИНС и пЭМГ-сигналом (а); 
соотношение между выходом ИНС и углом движения локтя (6) 


Угловое положение локтя оценивалось второй ИНС, качество распознавания которой варьировалось в зави- 
симости от набора входных параметров и количества нейронов скрытого слоя. Результаты моделирования процесса 
распознавания для второй ИНС приведены в табл. 3. 


Таблица 3 
Значение СКО распознавания углового положения локтя 





Кол-во нейронов | МЕ$Р,| 7С, | МЕЗР, МЕ$ЗР, |МЕБЗР, [ше,| МЕБЗР, ше, ГС, | МЕЗР, ие, ГС, 
в скрытом слое | Ге |1ое)еЁ те, ХС |Пие, 10зЛрей СС, 10ог)е! | 1овреь, КМ5 |10зре!, КМ$, 5ТР 





Среднеквадратическое отклонение 














5 1,610 | 1,791 | 1,589 1,567 1,630 1,528 1,545 
10 1,601 | 1,733 | 1,562 1,547 1,499 1,496 1,474 
15 1,586 | 1,744 | 1,552 1,549 1,571 1,556 1,511 
20 1,554 | 1,715 | 1,523 1,518 1,526 1,496 1,503 
































Выходной сигнал ИНС изменялся от 0 до 90, что соответствует углу поворота локтя 0-90 градусов. Наимень- 
шая ошибка распознавания получена при 10 нейронах скрытого слоя и использовании МЕБР, Пие ‚ СС, ог)ер, КМ5 во 
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входном векторе. Величина ошибки 1, 496 (2 %) может считаться малой. 

С помощью третьей ИНС определялся характер движения: +1 — сгибание локтя, 0 — отсутствие движения 
(пауза); —1 — разгибание локтя. На вход третьей МНС подавались сигналы: МЕ$Р, Пие, СС, орет, КМ, 5ТО. Резуль- 
таты моделирования работы (распознавания) ИНС приведены в табл. 4. 
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Значение СКО определения характера движения 





Таблица 4 




















ь МЕЗР, | 2С, |МЕБР, | МЕЗР, |МЕЗР, [ше, | МЕУЗР, Пие, СС, | МЕЗР, Пе, ГС, 
ОВО НоИрОНОЕ те |108)е! Шие, ХС ие, [ов)ей 7С, [0о—)е! | 10огре, КМ5 |1ов)ер КМ5, 5ТЬ 
в скрытом слое 
Среднеквадратическое отклонение 

э 0,2980 | 0,3235 | 0,2984 | 0,2992 0,2887 0,3593 0,3864 

10 0,2943 | 0,3742 | 0,3255 | 0,3298 0,2485 0,2295 0,2296 

15 0,2955 | 0,4491 | 0,3286 | 0,3290 0,3893 0,2987 0,2888 

20 0,2956 | 0,4413 | 0,3297 | 0,3328 0,3992 0,3385 0,3086 
































Несмотря на отсутствие груза в руках испытуемых, при сгибании — разгибании локтя формируются мышеч- 
ные усилия. Это объясняется тем, что необходимо перемещать массу самого предплечья и кисти. Поскольку мышеч- 
ное усилие создает вращающий момент, величина момента, необходимого для удержания предплечья в определенном 
положении, зависит от углового положения. Были рассмотрены зависимости развиваемой мышечной силы и момента 
от положения предплечья с учетом изменения размеров мышцы. 

Регистрация угла поворота, определение угловой скорости, длины мышцы бицепса и развиваемого момента в 
соответствии со схемой (см. рис. 1) выполнялись в 5йпийик системы МаНаБ. В псевдореальном масштабе времени 
рассчитывались: скорость движения, текущая длина мышцы и момент силы. 

Изменения длины бицепса и развиваемого момента при движении локтя без нагрузки в относительных еди- 
ницах приведены на рис. 3. Видно, что наибольшее сокращение мышцы соответствует максимальному моменту. 
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Рис. 3. Изменение длины бицепса и момента сил (мужчина, рост 178 см; 
объем бицепса 30/33 см; вес 82 кг). Движение локтя без нагрузки 
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Момент силы, развиваемый бицепсом (10), и его длина (11) определись в соответствии с рис. 1: 
В аа т 0(Е) 
— /азт0(0))2+(а+а соз 6)?” 


1 = \/(а$16(#))2 + (а +асоз0(2))?. 


(10) 
(р 








Очевидно, что пЭМГ-сигнал, фиксируемый при выполнении движений локтем, содержит информацию о па- 
раметрах движения и мышечных усилиях. Экспериментальные исследования с использованием обученных нейронных 
сетей позволили установить следующее. При определении угла поворота предплечья имеет место разброс (дисперсия) 
выходных значений нейронной сети (рис. 4). Для повышения качества распознавания после ИНС был использован 
фильтр Савицкого — Голея с порядком К = | и размером кадра } = 9. При этом результирующая величина СКО соста- 
вила 0,13, а график изменения выходного сигнала не имел разрывов первого рода. 
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Рис. 4. Угол поворота локтя (в относительных единицах). Показаны: выходной сигнал ИНС, 
сигнал датчика углового положения (гироскоп), результат фильтрации выходного сигнала ИНС 


Аналогичные экспериментальные исследования работы ИНС для распознавания характера движений разгиба- 
ния/сгибания позволили установить, что характер изменения выходного сигнала ИНС соответствует характеру изме- 
нения скорости движения предплечья (рис. 5). 

В силу нестационарности процессов, сопровождающих иннервацию мышечной активности и процессы в 
мышцах при выполнении движений, имеет место разброс выходных значений ИНС. Применение фильтра Савицкого 
— Голея (К = 1, }= 9) позволяет значительно улучшить качество выходного сигнала во всех фазах движения. Из полу- 
ченных результатов экспериментальных исследований и моделирования процессов распознавания можно сделать сле- 
дующее заключение. Применение обученной ИНС и последующая фильтрация ее выходного сигнала позволяют осу- 
ществлять оценку скорости выполнения движений локтем. 

1,0 


0,6 
0,4 
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Рис. 5. Распознавание характера движения предплечья. Показаны: выходной сигнал ИНС, 
результат его фильтрации, скорость предплечья (измерялась гироскопическим датчиком) 


Среднеквадратическая ошибка распознавания (оценка) угловой скорости составила 0,26, а после фильтрации 
— 0,20. 
Заключение. Выполненные экспериментальные исследования, моделирование и анализ полученных результатов по- 
казали, что наилучшие результаты распознавания наличия движения могут быть получены посредством ИНС, обу- 
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ченной с вектором входных параметров МЕЗР, [м!е. Также установлено, что с использованием ИНС и соответствую- 
щих процедур фильтрации на основе пЭМГ-сигналов мышечной активности возможно решение задач оценки пара- 
метров движения предплечья, в том числе для оценки угла поворота и оценки скорости движения. Среднеквадратиче- 
ская ошибка оценки угла и скорости в исследованиях не превышала 20 %. Таким образом, решены задачи, позволяю- 
щие разрабатывать системы биоэлектрического управления мехатронными устройствами. 
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